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During the 1996 multi-institutional Shelfbreak PRIMER experiment, low 
frequency sound sources were moored on the continental slope south of Cape Cod. These 
sources transmitted phase encoded tomography signals which were monitored by vertical- 
line hydrophone arrays moored on the continental shelf. The measured signals were 
processed for the acoustic modal arrivals and their variability in time. The processing 
entailed pulse compression, coherent averaging, local sound-speed profile updates and an 
application of the Chiu-Miller-Lynch model-based modal beamforming technique. In this 
thesis, the signal processing procedure is discussed and the modal arrival estimates are 
examined. The model-based estimates are found to be of high quality, with all 
propagating modes individually resolved. This unambiguous separation of the high 
modes cannot be achieved using simple least-squares techniques because of under 
sampling. The temporal variability of the modal amplitudes and travel times are found to 
be related to ocean processes that are unique to the shelf-slope littoral environment. 
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I. INTRODUCTION 



A. THE SHELFBREAK PRIMER EXPERIMENT 

The Shelfbreak PRIMER field study was a multi-institutional, joint effort by the 
Naval Postgraduate School (NPS), the Woods Hole Oceanographic Institution (WHOI), 
the University of Rhode Island, and Harvard University. As the oceanographic conditions 
at the shelfbreak differ drastically between seasons, the field work included two intensive 
three-week experiments, one in July 1996 (summer) and the other in February 1997 
(winter). The acoustic goal was to understand the propagation of sound from the 
continental slope to the continental shelf, including the effects of shelfbreak frontal 
features. The oceanographic goal was to characterize the mean and seasonal changes of 
the ocean circulation and determine the kinematics and dynamics of the shelfbreak frontal 
processes. [Refs. 1-2] 

The Shelfbreak PRIMER study area is shown in Figure 1. The area is south of 
Cape Cod in the vicinity of the continental shelfbreak. The experimental location 




Figure 1 . Geographic location of the 1996 Shelfbreak PRIMER experiment. 
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provided a complex coastal environment for studying sound propagation, from the 
continental slope to the continental shelf which is influenced by oceanographic and 
geological processes related to the shelfbreak. The shelf waters are cold and fresh, and 
vary considerably with season and location, while the warmer slope water is influenced 
by the Gulf Stream. The study area is marked by a significant ocean front at the shelf- 
break, the presence of strong internal tides, solitary waves (solitons), and complex flows, 
eddies and meanders [Refs. 1-4]. 

The Shelfbreak PRIMER field study was a cross-discipline physical oceanography 
and underwater acoustics experiment. In an effort to provide as much environmental data 
as possible a large suite of oceanographic and acoustic instruments were deployed. 

Figure 2 shows the principal acoustic moorings deployed during the summer experiment. 




40 -72 



Figure 2. The mooring configuration of the Shelf-break PRIMER summer 
experiment. The VLA’s (VLA1 and VLA2) were moored on the continental 
shelf, while the sound sources (SSI, SS2, and SS3) were moored on the slope. 
The highly variable bathymetry of the study area is also displayed. 
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The primary sensors for the acoustic component of the experiment were two 
vertical-line hydrophone arrays (VLA’s). Each array consisted of 16 hydrophone 
elements that spanned the lower two-thirds of the water column. The emitters used for 
the experiment were one 224 Hz tomography source with a bandwidth of 32 Hz, and two 
400 Hz tomography transceivers, each with a bandwidth of 100 Hz. This thesis study 
focused on the modal propagation along the western acoustic track. The transmitted 
signals from the 224 Hz (SSI), and 400 Hz (SS2) sources, co-located at the South-West 
comer, and received by the North-West vertical line array (VLA1) were processed and 
analyzed within the framework of acoustic normal mode theory. [Refs. 5-8] 



B. THESIS OBJECTIVES, APPROACH, AND OUTLINE 

The objective of this thesis is to (1) optimally process the acoustic VLA data for 
modal arrivals (magnitudes and travel times), and (2) relate the observed variability in the 
modal arrivals to oceanographic processes. The signal processing can be subdivided into 
pulse compression and coherent averaging (temporal processing), and modal 
beamforming (space-time processing). The transmitted signals were multiple periods of 
phase-modulated pseudo-random digits. A coherency study was performed to determine 
the optimal averaging time and optimum averaging gain. Once the individual 
hydrophone pulse arrival structures were processed, a model-based modal beamformer 
was applied to the data [Ref. 9], using updated sound-speed profiles local to the VLA. 
The variability analysis entailed examining the time series modal magnitude and travel 
time extracted from the modal-beamformed data. 

The remainder of this thesis consists of three chapters. Chapter II describes the 
signal processing methods used to process the VLA hydrophone data for the modal 
arrival structure. Chapter III presents and discusses the data processing results, and 
discuses the observed acoustic variability in relation to the shelfbreak oceanographic 
processes. Chapter IV presents the conclusions of this study. 
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II. SIGNAL PROCESSING 



A. PULSE COMPRESSION AND COHERENT AVERAGING 

1. Source Signal Characteristics 

The signal design and processing of ocean acoustic tomography has traditionally 
used the approach pioneered by T. Birdsall and K. Metzger. The Birdsall-Metzger 
method uses phase encoded periodic signals to eliminate self-clutter. These signals are 
easily generated and processed using digital techniques [Ref. 10, 1 1]. By modulating the 
signal’s phase with a pseudo-random binary sequence, a deterministic narrow pulse can 
be realized upon matched filtering. The sequence can also be transmitted repeatedly to 
achieve additional signal processing gain at the receiver by averaging consecutive 
sequences. 

Maximal-length binary sequences (m-sequences) have been used to generate a 
pseudo-random signal that can be decoded with no side-lobes in the time domain, 
provided that the sequence period is much greater than the arrival spread of the received 
signal [Ref. 12]. Figure 3 shows how the m-sequences are generated using linear shift 
register techniques. The shift register taps identify the sequence’s “law” value, typically 




Figure 3. Shift register generation of the pseudo-random M-sequence with a law 
value of 103. This m-sequence was used in the 224 Hz sound source (SSI). 
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represented as an octal value, which defines the binary code [Ref. 13]. For the 224 Hz 
tomography source, SSI, the m-sequence law used was 103 8 , or 100001 1 2 . The exact m- 
sequence generated for the experiment used an initial shift-register load of one. The same 
m-sequence (except for a time delay associate with a random seed of the shift register) 
would be generated with any random initial load, since an m-sequence is only determined 
by its law value. The source signal is generated by phase encoding the binary sequence to 
a suitable carrier for acoustic transmission, with each m-sequence digit encoded by a 
fixed number of carrier cycles. It is the number of carrier cycles per m-sequence digit 
that defines the temporal extent of the pulse of the output of a matched filter, which is the 
inverse of the modulating filter. 



. . . . 


SSI 


SS2 


SS3 


Center Frequency, f 0 


224 Hz 


400 Hz 


Full Bandwidth 


32 Hz 


100 Hz 


m-sequence law (octal) 


103 


1533 


1175 


m-sequence length, L 


63 digits 


5 1 1 digits 


digit width (resolution) 


14 cycles = 62.5 msec 


4 cycles =10 msec 


m-sequence period, T s 


3.9375 sec 


5.1 1 sec 


m-sequences transmitted 


30 sequences 


48 sequences 


transmission length 


118.125 sec 


245.28 sec 



Table 1. Source Signal characteristics for pulse compression in Shelf-break PRIMER. 

2. M-Sequence Period Mismatch 

In most previous tomography experiments, the m-sequence transmissions were 
typically periodically averaged to reduce both the size of the data and to attain coherent 
averaging gain. The Shelfbreak PRIMER VLA1 sample rate was hardware configured to 
be 1395.0892 Hz (divided down from a 5 MHZ reference frequency). The non-integer 
sample rate introduced minor difficulties in the periodic averaging of the received signal 
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and in the construction of the matched filter. These difficulties as well as their solutions 
are discussed next. 

The m-sequence period, T s , of each transmitted signal is listed in Table 1. Since 
the VLA digitized the received signal at a non-integer sample rate, the period of the 
received m-sequence was always mismatched to that of the transmitted sequence due to 
the discrete nature of the digitization process. To illustrate this mismatch, let us use the 
224 Hz transmission as an example. At a sampling frequency of 1395.089 Hz, the length 
of the received sequence is truncated to 5493 samples, while the true period length of the 
transmitted sequence is 5493.16406 samples. This mismatch represents a time skew of 
1 17.6 ps per m-sequence. Since the source signal consists of thirty consecutive 
sequences in each transmission, a time-skew of 3.528 msec can be found over the 
transmission without compensating for this mismatch. For the 224 Hz signal, with a digit 
resolution of 62.5 ms, a skew of 3.5 ms is negligible. In contrast, the 400 Hz signal is 
much more sensitive to the mismatch, introducing a time skew much greater than the m- 
sequence digit width of 10 ms. This skew can cause severe degradation in the subsequent 
coherent averaging process without a proper correction. Errors introduced by the period 
mismatch without correction are summarized in Table 2. 





SSI 


SS2 


SS3 


m-sequence length 


5493.16406 samples 


7128.90625 samples 


perceived length by VLA 


5493 samples 


7128 samples 


single m-sequence skew 


117.6 ps/m-sequence 


649.6 ps/m-sequence 


skew over one transmission 


3.528 msec 


31.18 msec 


of multiple sequences 









Table 2. Summary of error due to uncorrected period mismatch. 
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Since the skew is based solely on sample rate selection and signal period, it can 
easily be corrected by applying a phase shift of 






fix(T s F) 

F 

s 



( 2 . 1 ) 



to each perceived period of the received sequences to re-align the received signal to that 
of the transmitted period. The ‘fix’ function used in Eq. 2.1 rounds the argument to the 
nearest integer, towards zero. 

Figure 4 shows both the uncorrected and corrected demodulated envelopes of a 
received 400 Hz tomography signal containing 48 m-sequence periods. The uncorrected 



uncorrecied 400 Hz tomography arrival structure 




27.65 27.7 27.75 27 8 27.85 27.9 27.95 28 28.05 28.1 20.15 

travel time (sec) 



phase corrected anrval structure 




27.65 27.7 27.75 27.8 27.85 27.9 27 95 28 28.05 28.1 28.15 

travel time (sec) _ 



Figure 4. Uncorrected (top) and corrected (bottom) 
period-mismatched arrival structure of 48 consecutive 
m-sequences in a single 400 Hz transmission for the 
1996 year-day 206.3231. 
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pulse arrival structure (top panel) shows a clear time skew that can severely degrade 
coherent averaging gain without compensation. 



3. Coherent Averaging 

The ocean environment is highly variable at the New England continental shelf, 
particularly in the vicinity of the shelfbreak. The ocean variability caused significant 
signal fluctuations during the course of the Shelfbreak PRIMER experiment; this 
variability has limited how many arrival structures could be coherently averaged to attain 
the maximum averaging gain. The time duration over which one can coherently average 
multiple arrival structures to increase the signal-to-noise ratio before signal saturation or 
degradation occurs can be thought of as the inherent coherency time of the ocean sound 
channel. 

A temporal coherency study was performed over a transmission series to 




coefficients of the 400 Hz pulse 
responses measured at different 
times with five reference 
responses. 
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determine the maximum number of arrival structures that could be coherently averaged 
before degradation in the pulse resolution was seen in the arrival structure. The 
coherency time can be estimated by cross-correlating the arrival structures, i.e., pulse 
responses, observed at different times with one at a reference time. Figure 5 shows the 
correlation results for both the top and bottom hydrophone arrivals for the 400 Hz signal. 
Different reference times were used to study the stationarity of the coherence. 

Defining the limiting correlation value, below which the signal decorrelates 
significantly, to be 0.5 a coherence time of approximately 60 seconds is observed. This 
coherence time is consistent with that estimated by Chiu, et. al. [Ref. 4] based on a 
computer modeling study accounting for the effects due to large scale amplitude internal 
waves which have been observed frequently in the same geographic area. Note that the 
coherence time is not exactly depth independent and stationary, but exhibits small 
variations. 

Further investigation was performed by analyzing the signal-to-noise ratio (SNR) 
of the cumulative coherent averages of a transmission. Figure 6 shows the cumulative 
coherently averaged arrival structures for a 400 Hz tomography signal. The signal power 
was calculated over the predicted signal multi-path spread which occupies the time 
interval from 27.56 to 28.7 seconds relative to the beginning of the signal’s transmission. 
The noise power was calculated over the tail of the arrival structure segment, between 
28.85 to 32.53 seconds. This SNR value was calculated for the first arrival structure, the 
coherent average of the first two arrival structures, and so on, up to the coherent average 
of 48 consecutive arrival structures. The theoretical coherent averaging gain in a static 
sound channel is 101og(N), where N is the number of pulse responses being averaged. It 
is expected that the experimental result should approximate this theoretical gain as long 
as the ocean sound channel remains coherent. Figure 6 shows the experimental SNR gain 
and how it compares to the theoretical gain. The measured SNR gain shows a saturation 
at approximately 50 seconds, which is consistent with the results of the cross-correlation 
analysis. Figure 7 shows the cumulative coherently averaged pulse arrival structures as a 
function of the averaging time. 
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Theoretical gain (SNR) for the 400 Hz 
cumulative coherent averaged pulse 
responses for the top and bottom 
hydrophones of the VLA. 



Cor»(*ftl Av«r*0*» t<Tp twnTT 




27 4 27 6 271 2a 282 214 28 6 28 8 



CumuiAl** Cotwrvnl A v«tgM lar phon* 16 




27 4 27.8 27 J 28 282 28 4 28 6 28 1 



Figure 7. Cumulative coherent 
averaged 400 Hz pulse arrival 
structures shown vs. the averaging 
time for the top and bottom VLA 
hydrophones. 



It is important to note that although the SNR gain saturates beyond 50 seconds (10 
pulse responses) of coherent averaging, the individual pulses actually degrade beyond 
that. This implies that the number of pulse responses that can be coherently averaged 
before loss of pulse integrity is limited to ten. 



B. MODEL-BASED MODAL BEAMFORMING 



The ability to separate the vertical sound-field structure into individual normal 
modes for the purpose of ocean acoustic tomography has been demonstrated using an 
iterative model -based modal beamformer [Ref. 9]. This same technique was used to 
perform modal beamforming on the m-sequence removed, coherently averaged pulse 
responses measured by the North-West vertical line array, VLA1. For completeness this 
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technique will be outlined next. 

Computed by solving the Helmholtz Equation, the orthogonal normal modes can 
be fitted to the measured vertical sound pressure field at the hydrophone array to obtain 
their coefficients. With a continuous line array, the orthogonality of the modes is 
guaranteed. In practice, deployable arrays can’t span the entire water column due to 
fishing activity, hardware and environmental restrictions. Another limiting factor is 
hydrophone spacing. Without densely populated hydrophone elements, one ends up with 
an under-sampled pressure field. This under-sampling destroys the exactness of the 
orthogonality, which can result in mode-filter leakage of energy from one mode to 
another. [Ref. 9] 

The complex envelope of sound pressure p(z k ,t m ) sampled by each hydrophone in a 
VLA is related to the baseband pressure spectrum P(z k ,fj ) via a Fourier transform: 



M - 1 

nzjj > E 

m= 0 

k= 1,2,3,...,£ 
fj- j 



M At 

M , Af - M 

2 2 2 
t=t n +mAt 

m U 



( 2 . 2 ) 



where z k is the depth of the P h hydrophone in the VLA, t m denotes the discrete times at 
which the arrival structure is sampled. At is the sampling interval and the fs are the 
baseband frequencies with/ 0 =0 Hz. The interval MAr defines an M-point time segment 
stretching the duration of the arrival structure that begins after time t 0 According to full- 
wave theory, these baseband spectra, P(z k /j ), can be expressed as a linear combination of 
normal modes sampled at the different hydrophone depths. Accounting for the 
contamination by additive noise and denoting the sound pressure spectrum of the 
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complex envelope of noise at each phone by N(z k /j), the measured spectra can be 
expressed as 



p(z„, /p=E wyv fsF c )+N(z t , /p . 

n k=l, 2,3,... ,K 

. M , M _ M 

2 2 2 



(2-3) 



where N is the number of modes to be used in the beamforming and F c is the carrier or 
center frequency of the transmitted bandpass signal. Weighted by their frequency 
dependent coefficients A n in the sum, the local normal modes Z n are functions of both 
depth, z, and frequency, f=f c +f P that satisfy the following governing equation and 
boundary conditions: 



dz 2 P(z) dz dz c{z) 

Z n { 0, f)= 0 

Z n {h',f)=Z n {h;,f), 1 = 1,2,...,/ 

1 d z n (h-,j)=—L- A z n {h;,f), 1=1,2,...,/ 



P(h.) dz 

4 f>=°- 

dz 



P iK ) 



dz 



(2.4) 



Additionally, the mode functions at each frequency possess the property of orthogonality: 



/^) Z,(z ’ z " fe- - 0 * ' 6 - ■ < 2 - 5) 



In Eq. 2.4, k n is the horizontal wave number (or eigenvalue) of the nth mode, h is the total 
thickness of the acoustic waveguide that contains a total number of internal interfaces, h i 
is the depth of the 1 th interface (with i=l denoting the sea floor and i>l denoting the 
sediment-sediment interfaces below), and superscripts + and' denote depths just below 
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and above h p respectively. Given the sound-speed profile c(z), density profile p(z), and 
the depths of the water-sediment and sediment-sediment interfaces at the VLA location, 
the local normal modes and their eigenvalues can be determined by solving Eq. 2.4 and 
normalized using Eq 2.5. [Ref. 9, 14] 

Once the normal mode functions are known and the time series measured by the 
VLA have been converted to baseband spectra, the major task in model based modal 
beamforming is to estimate the complex modal amplitudes A n (i.e., the weighting 
coefficients) in the presence of “noise.” Depending on the number of modes that are 
processed for, the meaning of noise shall be expanded to include model error, i.e., the 
signals in the higher modes that are not included. Once an estimate A n (fj) of A n (fj) is 
obtained, an estimate a n (tj of the modal arrival structure a n (tj in time can be constructed 
using an inverse Fourier transform 

“„«J * E A.V? 

. M , 

J "~ ( 2 . 6 ) 

m= 0,1,2,. ..,M-1, 

A/=— 

M At 



To ensure that the mode beamformer is computationally tractable, even for large 
bandwidth signals, Chiu et al. [Ref. 9] formulated a procedure that operates on the 
baseband spectral data one frequency at a time. The measurement equation (Eq. 2.3) at 
each baseband frequency,^, can be rewritten in matrix- vector notation as 

Pj = ZjAj+ Nj ( 2 . 7 ) 

where the subscript is used to denote quantities associated with the baseband 
frequency./-. Associating Eq. 2.7 with Eq. 2.3; Pj is identified as a column vector of 
length K, containing the measured spectral values of the complex envelopes of sound 
pressure at each of the hydrophone depths, Z } is a K x N “mode shape” matrix, whose 
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columns are the “mode vectors,” each containing the depth samples of a specific mode 
function evaluated at the hydrophone depths, Ay is a column vector of length N denoting 
the unknown complex modal amplitudes A n at frequency j, and Nj is a column vector of 
length K representing additive noise. 

The “best” estimate Ay of Ay may be defined as the one that has the least mean 
square error and is linear with the data Pj , i.e.. 



Ay = OjPj. 



( 2 . 8 ) 



Under the assumption that second moment statistical models for A y and Pj are 
available, the Gauss-Markov theorem asserts that such an estimate is unique and that the 
corresponding estimator, for linear observations, is given by 

Oj = <A ^ T > Zj T (Zj <A jAj T > Z. T + <NjNj r >) “ 1 
or (2.9) 

Oj = (<AyAy V 1 +Z j T <N j N j T >' Zj)-' . 

where the brackets <•••> denote the expected value or first moment statistics. 

The corresponding error covariance matrix of the solution can be expressed as 

<(Ay-Ay)(Ay-Ay) r > = <A jAj T > - Of <A /l j T > Zy T f . (2.10) [Ref. 9] 

The Gauss-Markov estimator allows for an a priori solution covariance, <AjA?>. 
It is through this covariance that one can inject additional constraints, obtained by 
propagation modeling, using basic environmental information, to improve the 
beamformer results. For example, with sufficient environmental information (i.e., sound 
speed, bottom density, surface roughness, etc.) a propagation model can estimate the 
received modal amplitudes and the time spread between the propagating modes. Given 
an adequate model prediction A n (fj) , of the complex modal amplitudes A n (fj), we can 
express their relationship by 
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( 2 . 11 ) 



wwy 2 **' a m- 

where G n (fj ) corresponds to a distortion and 6t n a translation [Ref. 9]. 

We can calculate a model for the solution covariance (Eq. 2.14) if we apply the 
following assumptions: 

1. Translation and distortion of the predicted modal arrival structure are 
statistically independent. This is true if the dominant ocean processes that 
cause translation are different from those causing the distortion (i.e., 
processes of different time scales). 

2. Both distortion and translation of the arrival structure of different modes 
are statistically independent. This assumption is, in general, inexact 
because different modes do not traverse non-overlapping sections of the 
water column, but is not a poor assumption when exact knowledge of such 
correlations are absent. 

3. The second moments of the modal distortion factors are given by 

<\G n (fp\ 2> = (1 +e n ) 2 , (2.12) 

where the e n s are the modal-energy mismatch parameters. When £„=0, the 
ensemble average of the energy levels of the actual modal arrivals is 
identical to the predicted energy level. As e n deviates from zero, the 
predicted and the ensemble averaged levels become dissimilar. 

4. The time translations have Gaussian statistics such that their probability 
density functions are given by 

P(K>~ — e 2 "'' • (2I3) 

\] 2no t 
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where o is the standard deviation of the arrival time of mode n from that 

1 n 

predicted by the model. [Ref. 9] 






(1 ♦eX(f>4.V / >. 



m=n 



^2 ^ 

exp[-(27t^.) 2 — - }Tjfj)A’(fj), m*n 



(2.14) 



This analytic covariance model contains statistical parameters for mode arrival 
time changes, and mode amplitude variations. The refinement of these statistical 
parameters along with the updated predictions of A n (from previous iterations) is key to 
achieving improved mode beamformer performance. If the complex modal amplitudes 
predicted by the model are explicitly written out (Eq. 2. 15) in terms of amplitude and 

arrival time, t , we can gain some insight into the type of additional constraints that are 



being applied to this model-based modal beamformer. It follows from Eqs. 2.14 and 
2.15, that 



<A n (f j )A^{f j )>=e 



4 2 



(a, ta \ > 

} 1 n 1 m 






■av/i.ij 



m*n 



(2.16) 



Equation 2.16 reveals that the additional information contains the predicted relative 
arrival times between the different modes as well as the uncertainty of this 

information. The uncertainty placed on the relative arrival time information is imposed 
as the expected variance of the predicted arrival time. It is important to note that no 
constraints have been placed on the absolute arrival times. If constraints were placed on 
the absolute travel times that were too stringent, any error between the modeled anchor 
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positions and the actual experimental configuration would result in significant error in the 
modeformed arrival structures. [Ref. 9] 

C. SOUND SPEED PROFILE UPDATES 

The estimation of the modal amplitudes (i.e., modal coefficients) requires that the 
normal modes at the location of the VLA be known. These normal modes depend 
strongly on the local sound speed profile. Because of the complex ocean processes 
associated with the shelfbreak, such as frontal meandering, eddy activities, internal tides, 
and solitary internal waves, the sound speed profile and hence the normal modes 
themselves are extremely dynamic with strong spatial and temporal changes. One can no 
longer afford to use a nominal sound speed profile to compute the normal modes and 
assume that the normal modes are static throughout the experiment. Such an assumption 
can result in a wrong set of normal modes for a given time, which when used for the 
modal beamforming can cause large modal sidelobes, making it difficult to resolve the 
individual modes. 

In order to account for the dynamics of the normal modes at VLA1, we must 
obtain an accurate description of the sound speed profile as a function of both depth and 
time. This description was accomplished using the temperature time series obtained by 
four Brankner thermistor pods mounted at four depths, 24, 30, 50, and 70 m on VLA1, as 
well as salinity-temperature-depth data sampled by a surveying SeaSoar. A SeaSoar is a 
towed oceanographic sensor moving up and down the water column, with the horizontal 
resolution determined by the towing speed. The underlying logic was that the Brankner 
measurements contain information on die temporal evolution, while the SeaSoar data 
contain information on the vertical structure of the ocean’s temperature. 

SeaSoar data sampled near VLA1 were first extracted to construct the local 
vertical structure of the ocean temperature using an empirical orthogonal mode (EOF) 
analysis [Ref. 14]. The Brankner temperature time series were then projected onto the 
EOFs and their temporal variations are shown in Figure 8. Finally, the analyzed 
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temperature profiles were converted into sound speed profiles using a standard empirical 
formula to facilitate an update, i.e. a recomputation, of the normal modes at each 
transmission time. 




Figure 8. The temporal evolution and spatial changes of the temperature profile at the 
VLA1 site. This continuous ocean temperature field was analyzed from data sampled by 
four Brankner temperature pods and a surveying SeaSoar vehicle. 
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III. RESULTS AND DISCUSSION 



A. MODAL ARRIVAL ESTIMATES 

The Northwest vertical-line hydrophone array (VLA1) collected acoustic data 
between 24 July (yearday 206) and 2 August (yearday 215), 1996. During this time the 
Southwest 224 Hz tomography source (SSI) transmitted a phase-encoded signal every 
five minutes, while the Southwest 400 Hz source (SS2) transmitted every fifteen minutes. 
The sampling done by VLA1 has two time gaps during the entire experiment. These gaps 
in the data time series were caused by hardware problems; a blown fuse caused a twenty 
hours of data to be lost, and a disk drive problem caused another nine hours of data loss. 
Another 18-hour gap can be found in the modal beamforming results between 1996 
yearday 209.6 to 210.4. During this time period, the sampling rate of VLA1 was 
increased to accommodate the broadband, explosive charge experiment which was 
another component of the PRIMER experiment. This higher rate data has not been 
beamformed at this time. 

There are nine propagating modes at 224 Hz and 17 propagating modes at 400 Hz 
at the VLA1 location. The model-based mode beamformer was applied to the pulse 
arrival structures measured by the 16 omni-directional hydrophones of the VLA to 
estimate the variability of the amplitudes and travel times of these propagating modes. 

To demonstrate the improvement in quality of the model-based results, a simple least- 
squares solution [Refs. 15, 16] to Equation 2.7 was also computed to establish a baseline 
for comparison. The estimated modal arrivals constitute a fundamental data set for an 
inverse mapping of the coastal ocean temperature [Ref. 17]. The inversion of the travel 
times is known as tomography, a subject beyond the scope of this thesis research. 
Hereafter, the simple least-squares method will be referred as LS, whereas the model- 
based method will be referred as MB. The LS method can be interpreted as one without 
the time-separation and amplitude constraints provided initially by the model-based 
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procedure described in the previous chapter. 

Figure 9 shows the arrival structures of the lowest mode (i.e., Mode 1) estimated 
by MB for both the 224 Hz and 400 Hz signals. When compared to the LS estimates 
displayed in Figure 10, it is evident that the MB estimate is of higher quality, particularly 




Figure 9. The arrival structure of Mode 1 estimated by the model- 
based beamformer as a function of transmission time for both the 
224 Hz (left panel) and 400 Hz (right panel) signals. 



for the 400 Hz signal. The modal sidelobes in the LS estimate for the 400 Hz signal are 
significant, resulting in multiple peaks, thus making the identification of the modal travel 
time ambiguous. The arrival structure estimated by MB is unimodal for each 
transmission. 

The modal sidelobes are more severe for the higher order modes in the LS 
estimate. To illustrate this, we use the estimated arrival structure of Mode 16 for the 400 
Hz transmissions as an example. The right panel of Figure 1 1 shows the LS estimate of 
the Mode 16 arrival structure. It is obvious that this higher mode is unresolved by LS 
because of the contamination or ambiguity presented by the severe sidelobes. The 
presence of modal sidelobes is inherent to spatial under sampling, inadequate hydrophone 
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Figure 10. The arrival structure of Mode 1 estimated by simple 
least squares as a function of transmission time for both the 224 Hz 
(left panel) and 400 Hz (right panel) signals. 




Figure 11. The arrival structure of Mode 16 estimated by simple 
least squares (left panel) and model based modal beamforming 
(right panel) for the 400 Hz signal. 
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spacing and aperture coverage of the water column. Higher modes and higher frequency 
energy are more affected by these sidelobes due to their shorter vertical wavelength and 
larger vertical span. The reduction of these sidelobes requires additional information, 
which can come from physical constraints that are meaningful, reasonable and refineable. 
This is the basic philosophy behind the MB method developed by Chiu et al. [Ref. 9], 
with progressive refinements based on the actual data. The right panel of Figure 1 1 
shows the MB estimate of a clearly resolved Mode 16. 

B. ACOUSTIC VARIABILITY 

As seen in figures 9, 10, and 11, both the magnitude and travel time of the 
individual normal modes can vary significantly over time. This acoustic variability 
consists of many time scales, from weeks to minutes, and is mode-number dependent. 

Due to the high-resolution oceanographic data available during the experimental period, a 
fairly good estimate of the environmental influences that caused this variability can be 
developed. 

On the close-to-weekly time scale, one can easily observe significant magnitude 
and travel time changes in the modal arrivals over all the modes, though they are much 
more apparent in the lower modes [Fig. 9]. The magnitude and travel time changes are 
well correlated; as the travel time increases, so does the magnitude of the pulse arrival. 

An increase in the travel time corresponds to an increase in the temperature along the 
transmission path. The key to understanding these results lies in the following questions: 
where along the path was the temperature increase, what caused the temperature to 
change, and why was the signal level affected? The SeaSoar, moored thermistors and 
infra-red satellite imagery provided the data needed to answer these questions. These 
oceanographic data show that there was a warm water filament that intruded over the 
transmission path twice over the course of the experiment. This warm filament was the 
remnant of a dissipated Gulf Stream, warm core ring. The core of the warm intrusions 
were isolated to the lower portion of the water column, and thus increased the sound 
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speed near the bottom. This sound speed increase resulted in a more upward-refracting 
sound channel, which reduced the signal attenuation due to bottom interactions (sediment 
absorption). After yearday 213, the warm filament dissipated and the lower water column 
on the shelf remained cold, and a much lower signal level is observed as a result of the 
much stronger interactions with the bottom. In general, the lower mode signal responses 
tend to be stronger because of the lesser bottom interactions. 

The arrival times of the individual modes and their variations over transmission 
time for both the 224 and 400 Hz signals were extracted by picking the peaks of the MB 
estimates. The resulting time series of modal travel times are shown in Figure 12. Low 
pass filtering was applied to remove random fluctuations due to internal waves with 
periods shorter than three hours. The observed travel times show that the low modes 
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Figure 12. Measured variability of the normal mode travel times 
over the duration of the summer experiment. There are nine time 
series plotted in the top panel, each corresponds to one of the nine 
propagating modes at 224 Hz. There are 16 time series plotted in 
the bottom panel, each corresponds to a propagating mode at 400 
Hz. The lower modes arrived sooner whereas the higher modes 
arrived later for both frequencies. 
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travel faster, and the higher modes travel slower. This result is expected in a shallow- 
water environment where the vertical extent of a mode increases with mode number. 
Higher modes are analogous to rays with steeper angles which travel through longer 
“raypaths” to reach the receiver. 

In addition to the weekly variability introduced by the intrusion of the warm 
filament, other significant acoustic fluctuations with periods less than a day can also be 
observed [Fig. 12]. Fluctuations due to internal, diurnal tides are dominant in the 
intermediate modes for the 400 Hz and in the higher modes at 224 Hz. These modes have 
a large expression in the main thermocline where the internal tidal energy is mostly 
confined. At shorter time scales of several hours, the acoustic fluctuations can be 
attributed to longer-period internal waves as well as the large-amplitude internal soliton 
wavepackets [Ref. 3]. 
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IV. CONCLUSIONS 



A. SUMMARY OF FINDINGS 

In this study, the model-based modal beamforming technique was applied to data 
recorded by a vertical line array (VLA) to estimate the arrival structure of the individual 
normal modes that propagated from the continental slope, through a shelfbreak front, to 
the continental shelf. The model-based estimate, when compared to a standard least- 
squares estimate, shows significant improvement. The improvements in the model-based 
estimate, in terms of filter resolution or sidelobes, are clearly seen for the higher modes. 
These modes are largely unresolved in the standard least-square estimate, while they are 
unambiguously resolved and identified in the model-based estimate. This results in a 
high quality data set of modal travel times to facilitate subsequent tomographic inverse 
analysis [Ref. 17]. 

Incorporation of the temporal changes of the normal modes was also found to be 
extremely important to maximizing the quality of the mode beamformer. This was 
accomplished by updating the sound speed profile based on an optimal interpolation 
analysis of the in-situ temperatures measured in the vicinity of the VLA. This update was 
then followed by re-computing the normal modes for each transmission time. 

An empirical study of the coherence time of this shelf-slope sound channel at the 
Shelfbreak PRIMER transmission frequencies was carried out to investigate how long the 
signal can be coherently averaged before saturation or degradation occurs. Based on data 
obtained by the top and bottom hydrophones of the VLA, this coherence time was 
estimated to be approximately 50 s, although mild nonstationarity and inhomogeneity 
were observed. This coherence time was used to determine the number of consecutive 
pulse responses that could be coherently averaged to maximize the SNR without 
spreading the response structure. 

A first look at the variability of the acoustic normal modes was taken using the 
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estimated modal arrival structure in conjunction with a high-resolution oceanographic 
data set. Significant fluctuations in modal travel times and magnitudes were found in 
multiple time scales due to oceanographic processes of different types such as warm-ring 
filament intrusions, internal tides, and internal waves. An in-depth study of the dynamics 
of these ocean processes and a quantification of their influences on the acoustic 
transmissions are currently going on. 

B. RECOMMENDATIONS FOR IMPROVEMENT 

Several improvements of the signal processing of the Shelfbreak PRIMER 
acoustic data are recommended. First, the bandwidth of the 400 Hz signal was 100 Hz. 
However, because of the model-based technique was developed only for narrow-band 
signals, only a sub-band of 32 Hz was processed. There are two options to proceed. One 
option is to modify the modal beamforming method to accommodate for the 100 Hz 
bandwidth. This can be accomplished by dividing the 100 Hz band into several narrower 
bands and treat each sub-band independently in the model-based mode beamformer. One 
can then combine the estimates for each non-overlapping sub-bands together in the 
inverse Fourier transform to construct the modal arrival structures of a higher resolution 
pulse. The other option is not to combine the sub-band estimates but treat them 
individually when constructing the time domain modal arrival structure. The former 
option has the advantage of giving more accurate travel time estimates, whereas the latter 
has the advantage of giving more travel time data associated with modes at different 
frequency bands. Modes at different frequencies sample the ocean slightly differently and 
may help to improve the resolution maps of the ocean frontal structure in a tomographic 
inverse analysis. The two options need to be evaluated and compared. 

By examining the acoustic data from the individual hydrophones of the VLA, the 
signal phases were observed to roll linearly or quadratically over each transmission. The 
phase rolls were observed to be stronger for the 400 Hz signals than for the 224 Hz 
signals and is consistent with the hypothesis of a Doppler shift induced by ocean currents. 
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It is therefore useful to examine if the Doppler can be compensated to improve the 
coherence time of the sound channel. 

A simulation study of the effect of the VLA tilt on the performance of the mode 
beamformer is also warranted. A large syntactic foam subsurface buoy with 2,000 lb 
buoyancy was used to moor the VLA with an expected tilt of less than one degree. It is 
useful to examine the degradation of the mode beamformer performance as a function of 
frequency. Although the results at 224 and 400 Hz may not be significantly affected, 
future shallow water experiments that will use higher frequencies can benefit from this 
study. An outstanding question is: At what frequency must the VLA tilt be monitored 
and what is the desired accuracy? The answer to this will be useful for design purposes. 

Improvements in the processed results could have been achieved if the tilt sensors 
on the VLA had functioned as expected. At 400 Hz, an array tilt of 1 .0 will result in a 
quarter wavelength phase shift in the data across the array aperture of 52.5 meters, 
severely affecting the results. Statistical estimation was performed on this data set to 
correct for the lack of tilt information resulting in additional uncertainty variances on the 
final results. 

Another low-cost, simple addition to the experiment configuration would be to 
increase the number of temperature pods used along the vertical array. Having the 
increased environmental information along the highly variable shelf-break allowed us to 
account for the temporal variability when calculating the vertical mode structures. 
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